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ABSTRACT 

Effective temperatures and radii for 92 planet-liosting stars as determined 
from tlie InfraRed Flux Method (IRFM) are presented and compared with those 
given by other authors using different approaches. The IRFM temperatures we 
have derived are systematically lower than those determined from the spectro- 
scopic condition of excitation equilibrium, the mean difference being as large 
as 110 K. They are, however, consistent with previous IRFM studies and with 
the colors derived from Kurucz and MARCS model atmospheres. Comparison 
with direct measurements of stellar diameters for 7 dwarf stars, which approxi- 
mately cover the range of temperatures of the planet-hosting stars, suggest that 
the IRFM radii and temperatures are reliable in an absolute scale. A better 
understanding of the fundamental properties of the stars with planets will be 
achieved once this discrepancy between the IRFM and the spectroscopic temper- 
ature scales is resolved. 

Subject headings: stars: fundamental parameters — stars: planetary systems 
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1. Introduction 

Accurate determination of the fundamental stellar parameters for the planetary host 
stars is required to improve our knowledge of the physical properties of the extrasolar planets. 
The temperatures and luminosities of these stars are used to obtain their masses and radii, 
which are the stellar parameters required to determine the masses and radii of the planets 
from radial velocity and transit observations, respectively. 

The temperatures of the planetary host stars have been spectroscopically determined 
from the excitation equihbrium of Fe I lines (e.g. Heiter & Luck 2003; Gonzalez et al. 2001; 
Laws et al. 2003; Santos, Israehan, & Mayor 2004). The most complete study of this kind 
is that of Santos et al. (2004, SIM04), who determined temperatures of 98 of them. Using 
a different approach, Ribas et al. (2003, R03) determined the temperatures of these stars 
from infrared photometry. Both results are similar and, in fact, consistent with most of the 
previous studies, with the exception of the Fe I temperatures of Sadakane et al. (2002) and 
Takeda et al. (2002), which are about 100 K lower. 

To date, it has not been verified whether the spectroscopic temperatures are consistent 
with those determined from stellar angular diameters. This is basically due to the lack of 
dwarf stars with measured angular diameters. However, recent interferometric measurements 
(Kervella et al. 2003, 2004; Lane, Boden, & Kulkarni 2001; Pijpers et al. 2003) and transit 
observations (Brown et al. 2001) allow one to test this now. Here we show that our IRFM 
temperatures are consistent with direct measurements of angular diameters. Then we apply 
the IRFM to obtain the temperatures of 92 of the known planet-hosting stars and compare 
them with those obtained by other methods. 

2. Determination of temperatures and radii 

Our work is based on the InfraRed Flux Method, which compares the quotient between 
the observed bolometric flux and the flux at a given wavelength in the infrared with the 
quotient between crTjj and the emergent flux in the infrared predicted by models to get the 
effective temperature and angular diameter of a star (Blackwell, Petford, & ShaUis 1980). 

The stellar atmospheric parameters required for our implementation of the IRFM have 
been obtained from an updated version of the Cayrel de Strobel, Soubiran, & Ralite (2001) 
catalogue and from SIM04. The good agreement between common [Fe/H] values in these 
two lists allowed us to use them together. The bolometric fluxes were calculated from the 
cahbration of Alonso, Arribas, & Martmez-Roger (1995, AAM95) and the photometry was 
taken mainly from 2MASS. Spectral energy distributions from Kurucz models were adopted. 
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The 2MASS photometry needed to be transformed into TCS colors given the differ- 
ences found between the various infrared systems (Carpenter 2001). From a sample of main 
sequence stars with TCS photometry in Alonso, Arribas, & Martinez-Roger (1994b) and 
reliable 2MASS photometry we found the following relations: 



with standard deviations of a = 0.038 (A^ = 104), a = 0.039 (A^ = 103) and a = 0.035 
(N — 107), respectively (A" is the number of stars in every fit). The brightest planet- hosting 
stars are not included in our study given the large errors associated to their 2MASS pho- 
tometry. However, we were able to calculate the temperatures of 5 bright stars (HD 009826, 
HD 022049, HD 027442, HD 117176, and HD 137759) from the Johnson photometry avail- 
able in the General Catalogue of Photometric Data (Mermilliod, Mermilliod, & Hauck 1997, 



Apphcation of the IRFM results in three temperatures, one for each band: Tj, Th and 

Tk- In general, they were weighted averaged with the inverse of their errors to get the final 
Teff, that is, Teff = (J2iTi/^Ti)/(J2i '^/^Ti); with i = J,H,K. The internal error in TeS was 
then calculated from AT^s = 3/(^ - 1/ATj). Further details on our IRFM implementation 
will be given in a forthcoming paper (Ramirez & Melendez, in preparation). 

One of the critical ingredients of the IRFM is the absolute infrared flux calibration 
of the standard star. We have adopted that given by Alonso, Arribas, & Martfncz-Roger 
(1994a), which was derived by demanding that the IRFM temperatures (angular diameters) 
be well scaled to the direct (measured) ones, at least for giants. This seems to be also true 
for dwarfs, as it is shown in Fig. 1 and Table 1, where the angular diameters of 7 dwarf stars 
derived in our work are compared to those measured by Kervella et al. (2003, 2004), Lane 
et al. (2001), and Pijpers et al. (2003) from interferometric observations and by Brown et al. 
(2001) from transit observations. 

The temperatures for these 7 dwarfs were obtained in the same IRFM implementation. 
However, with the exception of HD 209458, the photometry was adopted from sources other 
than 2MASS, given that they are bright stars. For r Cet and GJ 105, TCS photometry was 
available (Alonso et al. 1994a). For Procyon and GJ380 we took the mean of the Johnson 
photometry available in the GCPD. Finally, for a Cen A and B, Johnson infrared magnitudes 
as observed by Engels et al. (1981) were used. Johnson photometry was transformed into 
the TCS system by using the transformation equations of Alonso et al. (1994a). 



Jtcs = ^2MASS + 0.001 - 0.049( J 
Htcs = ^^2MASS - 0.018 + 0.003( J 
i^TCS = i^2MASS - 0.014 + 0.034( J 



-f^)2MASS 
-f^)2MASS 
-f^)2MASS 



(la) 
(lb) 
(Ic) 



GCPD). 



Very large systematic errors in the adopted /boi scale may lead to a wrong T^s scale. Our 
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work adopts the fluxes as derived from the AAM95 cahbration, which is based on integrated 
UBVRIJHK photometry and a theoretical correction for the flux outside the U-K wavelength 
range. AAM95 showed their fluxes to be consistent with previous works as that of Blackwell 
& Petford (1991), which was based on spectrophotometric data, at least for solar-metallicity 
stars. They also showed that the observed differences were obviously due to metallicity 
effects. 

A comparison with the bolometric fluxes that are derived from the theoretical BC scale 
of Bessell, Castelli, & Plez (1998) for the stars in Table 1 shows a tendency with temperature. 
The fluxes from the theoretical BC's seem to be higher by about 3% at 4750 K but lower by 
about 1.5% at 6500 K (open circles in Fig. 2). Also, when comparing measured fluxes for 
35 dwarf stars from the Blackwell & Lynas-Gray (1998) work with the measured fluxes used 
by AAM95 to derive their cahbration, the differences show a similar behavior, in the sense 
that they increase with temperature but only from —1% to +1%. The flUed circles and error 
bars in Fig. 2 correspond to the mean difference and ± 1 cr bars in 500 K bins, respectively. 
These small discrepancies, however, do not affect our flnal results substantially as will be 
explained in Sect. 3. 

Measurements of /boi are only available for Procyon (Code et al. 1976, Beeckmans 
1977, Blackwell & Shaflis 1977, Smalley & Dworetsky 1995, Mozurkewich et al. 2003), r Get 
(Blackwell & Lynas-Gray 1994, 1998) and GJ 105 A (Blackwell & Lynas-Gray 1998). By 
comparing the mean values of these measurements with the /boi values we adopted, a mean 
difference of (0.20 ±1.73)% is found (the major part of the dispersion comes from GJ 105 A). 
If we include the fluxes measured by AAM95, the difference reduces to (0.08 ± 0.94)%. 
Since there is no apparent reason to exclude the AAM95 fluxes from the others, we may 
conclude that the adopted /boi scale in our work is consistent with the measurements and 
the dispersion is within the adopted mean error for the fluxes. The stars in Fig. 2 show 
the difference between the flux as derived from the AAM95 calibration and the mean of the 
measurements. 

The measured fluxes and the fluxes from the calibration of AAM95 are given in Table 1. 
They have been used to get the direct temperatures given in the same Table 1, the measured 
fluxes have been preferred, of course. For the 4 stars with no /boi's measTircd. we have used 
the temperatures from the calibration. This is a reasonable approximation given that the 
error in the flux is propagated to o-(Tcff) roughly as o"(/boi)/4 and that there is evidence for 
the calibration to be in good agreement with the measured fluxes. 

The stars in Table 1 cover the range of temperatures from 4000 K to 6600 K, which is 
approximately the range covered by the planet-hosting stars. The zero point of the IRFM 
Teff scale for the dwarf stars is in good agreement with the direct temperatures. Gonsidering 
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the best five measurements in Table 1, the mean difference {T^s™ — T^s) with 
a standard error of 24 K. Note that the angular diameters measured for GJ 105 A and 
HD 209458 have large errors and so they are not useful to constraint any T^s scale. In fact, 
within ±200 K, any existing T^e scale agree with these measurements. 

The temperatures derived in this work are given in Table 2. Also given in Table 2 are the 
radii of these stars in solar units, obtained from their Hipparcos parallaxes, the bolometric 
flux scale of AAM95 and the IRFM angular diameters. 

A Tes scale for the planet-hosting stars is obtained by using the {V — K) color index, 
where K is from 2MASS: 

= 0.460 + 0.313(y -K)- 0.025(V - Kf . (2) 

This formula was obtained from a fit to the data for the dwarf stars with reliable 2MASS 
photometry in the list. Its standard deviation is 43 K and its ranges of applicability are 
1.20 <{V -K) < 2.60, -0.50 < [Fe/H] < 0.45. A plot of the data and Eq. (2) is given in 
Fig. 3a. Since the metallicity range covered by these stars is short, the metallicity dependence 
of this relation could not be distinguished from the internal errors in the temperatures and 
so no additional terms containing [Fe/H] were included in the fit. 



3. Comparison with other studies 

Our calibration for the T^s vs {V — K) relation is in agreement with the colors derived 
from Kurucz models (M. Bcsscll, private communication), as it is shown in Fig. 3b, where 
the colors for [Fe/H] = —0.5 ,0.0 and +0.5 are plotted along with our calibration. The 
colors kindly provided by M. Bessell are in the Bcssell & Brett (1988) system and have 
been transformed into 2MASS colors by using a transformation equation given by Carpenter 
(2001). Although it is not shown in the figure, a good agreement is also found with the 
solar metaUicity Teflf vs {V — K) relation of Houdaschelt, Bell, & Sweigart (2000), which is 
based on MARCS models. On the other hand, the temperatures given by R03 and SIM04 
lead to a T^q scale that has to be shifted by about 110 K to better agree with them (squares 
and open circles in Fig. 3a). Note that a similar difference {T^^'^ — T^^^ ~ +139 K) was 
found by S1M04. Fig. 4a and 4b show the differences between our Tcfj's and those given by 
SIM04 and R03, respectively. The calculated mean differences are —115+83 K (91 stars) for 
the points in Fig. 4a and —108 + 53 K (79 stars) for Fig. 4b. The label 'solar (yfj-values' in 
Fig. 4a emphasizes the fact that SIM04 used these transition probabilities in their analysis, 
while 'synthetic VJHK photometry' in Fig. 4b stands for the basic idea of the R03 method. 
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The adopted /boi scale in our work can not be the only cause of the difference, and, in 
fact, we believe it is not even the main cause of it. The dotted lines in Fig. 2 illustrate how 
an offset in the derived IRFM temperatures corresponds to a systematic error in the adopted 
bolometric fluxes. They have been derived by assuming that the /boi's are the only source 
of error in the implementation and so have produced temperatures systematically lower and 
higher by 30 K compared to the absolute scale. The quantities A/boi are then the differences 
between the adopted /boi's and the real ones. Even if the AAM95 calibration is inaccurate, 
its effect on the Teff's can not be larger than about 30 K, and so it is very unlikely for it to 
be the main cause of the 110 K difference. 

The temperatures obtained from photometric calibrations based on previous IRFM stud- 
ies (e.g. Blackwell & Lynas-Gray 1998; Melendez & Ramirez 2003) agree very well with the 
present results. As an illustration, the IRFM Tes calibrations for the Geneva colors {B2 — G), 
{B2 — Vi) and the t parameter from Melendez & Ramirez (2003, MR03) were apphed and 
averaged to get the Tefj's and compare them with the present results in Fig. 4c. Also shown 
in Fig. 4c is the difference between the Tcfj's obtained from the Blackwell & Lynas-Gray 
(1998, BLG98) IRFM %s calibration for {B2 - Vi) and our %s's. 

Not all the temperatures obtained from the excitation equilibrium of Fe I lines are in 
disagreement with the IRFM. The temperatures obtained by Sadakane et al. (2002, S02) 
and Takeda et al. (2002, T02), for instance, are in reasonable agreement with the IRFM 
although a slight tendency with temperature is found (Fig. 4d), the later maybe due to 
non-LTE effects. This may also be the cause of the small slope observed in Fig. 4a, which is 
roughly similar to that in Fig. 4d. The main difference with the work of SIM04 seems to be 
the use of transition probabilities measured in laboratory instead of solar gf -values. There 
is also a larger coverage in excitation potential of Fe I lines more sensitive to temperature, 
SIM04 use 3 lines with Xexc < 2 eV whilst S02 and T02 use 12. Only two of these hues 
are in the blue region of the spectrum, where the continuum can not be accurately defined. 
However, all the lines employed by S02 and T02, in addition to have reliable experimental 
transition probabilities, were carefully selected and so are hardly affected by blending. 



4. Conclusions 

The temperatures obtained in this work for the stars with planets differ from those 
given by most of the other groups by about 110 K. In particular, there is a discrepancy with 
the excitation equilibrium temperatures of Santos et al. (2004) and the method of infrared 
photometry of Ribas et al. (2003). 
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On the other hand, the Teff vs {V — K) relation from model atmospheres is closer to our 
Tcff scale and the angular diameters measured by interfcromctry and transit observations for 
7 dwarf stars are well scaled to those derived from the IRFM. The bolometric flux calibration 
adopted (AAM95) seems to be in excellent agreement with the measurements reported in 
the literature (±1%), which strengthen the present results. The comparison with direct 
measurements suggest that the zero point of the IRFM Tgff scale is within 50 K, considering 
twice the standard error. New measurements of angular diameters for stars in the range 
4000 K< Teff <6000 K are encouraged to better define it. 

The 110 K difference makes the stars bigger according to the IRFM, which implies that 
the planetary radii obtained from transit observations would be larger if our T^q scale is 
adopted instead of the spectroscopic one. 

Metallicity determinations for the stars with planets taking into account the IRFM Tgfj 
scale would be interesting. Note, however, that large samples of stars both with and without 
detected giant planets have been analyzed in a homogeneous way, so it is unlikely that this 
will affect the well stated metaUicity enhancement of the planet-hosting stars, although small 
differences may be found in an absolute scale. 

We thank M. Bessell for providing colors for the complete set of Kurucz models and 
I. Ivans for her comments and suggestions. Significant improvement to the original manuscript 
was possible due to a very constructive refereeing process. This publication makes use of 
data products from the Two Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and IPAC/Caltech, founded by NASA and NSF. 
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Table 1. Cool Dwarfs with Measured Angular Diameters 



Star name 


HD 


^LD (mas) 


Rcf.'^ 


/bol 




fhof 






ymFM (K)c 


6m.FM (mas) 


Procyon 


61421 


5.448 ± 0.052 


1 


1.837E- 


'5 


1.822E- 


-5 


6553 ± 40 


6591 ± 73 


5.405 ± 0.131 


r Get 


10700 


1.971 ± 0.050 


2 


1.163E- 


-6 


1.154E- 


-6 


5473 ± 72 


5372 ± 65 


2.047 ± 0.054 


GJ 380 


88230 


1.306 ± 0.040* 


3 


1.373E- 


-7 






3962 ± 63 


3950 ± 161 


1.301 ± 0.107 


GJ 105 A 


16160 


0.941 ± 0.070 


3 


1.687E- 


-7 


1.725E- 


-7 


4917 ± 185 


4714 ± 67 


1.013 ± 0.031 




209458 


0.226 ± 0.015 


4 


2.279E- 


-8 






6049 ± 202 


5993 ± 71 


0.230 ± 0.006 


a Gen A 


128620 


8.511 ± 0.020 


5 


2.677E- 


-5 






5771 ± 23 


5759 ± 70 


8.548 ± 0.225 


a Gen B 


128621 


6.001 ± 0.034 


5 


8.653E- 


-6 






5182 ± 24 


5201 ± 65 


5.957 ± 0.160 



^(1) Kervella et al. (2004); (2) Pijpers et al. (2003); (3)Lane et al. (2001); (4) Brown et al. (2001); (5) Kervella et al. (2003). 

''Bolometric fluxes from the Alonso et al. (1995) calibration. Units are erg cm~^ s~^. 

''Mean of the measured bolometric fluxes (see the text for references). Units are erg cm~^ s~^. 

'^Direct temperatures calculated by using the measured bolometric fluxes or the adopted flux calibration. An error of 1.5% 

was adopted for the /toi values. 

"Temperatures from the present IRFM implementation. 
'Estimated from the UD diameter with a correction of 3%. 
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2 4 6 8 10 

^LD (mas) 

Fig. 1. — Comparison between measured angular diameters (^ld) and diameters from the 
IRFM (^irpm). 
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Fig. 2. — Difference between the ffuxes from tfie Alonso et al. (f995, AAM95) calibration 
and tlie fluxes from: Bessell et al. (1998) bolometric correction scale for the stars in Table 1 
(open circles), Blackwell & Lynas-Gray (1998) measurements for 35 stars in 500 K bins (filled 
circles, mean differences and ± 1 cr error bars are shown), and other measurements (see the 
text for references) for 3 stars in Table 1 (stars) . The dotted lines illustrate how an offset of 
±30 K in the derived IRFM Tes scale may be due to wrongly adopted bolometric fluxes. 
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Table 2. Temperatures and Radii of the Planet-Hosting Stars from the IRFM 



HD 




(K) 




R/R 


'0 




HD 


Tcff 


(K) 


R/Rq 


HD/BD 


Teff (K) 




R/R 


'0 




000142 


6152 


± 


77 


1, 


.467 


± 


0. 


.032 


065216 


5561 


± 


66 


0.900 


± 0.036 


143761 


5754 


± 


70 


1, 


.341 


± 


0. 


.029 


001237 


5512 


± 


68 


0. 


.832 


± 


0. 


.028 


068988 


5778 


± 


70 


1.209 


± 0.062 


145675 


5129 


± 


67 


1. 


.022 




0. 


.030 


002039 


5847 


± 


69 


1, 


.251 


± 


0. 


.105 


070642 


5620 


± 


67 


1.023 


± 0.031 


147513 


5852 


± 


62 





.950 


± 


0. 


.026 


003651 


5264 


± 


67 


0, 


.844 


± 


0. 


.029 


072659 


5798 


± 


69 


1.481 


± 0.059 


150706 


5792 


± 


69 





.964 


± 


0. 


.030 


004203 


5546 


± 


68 


1, 


.418 




0. 


102 


073256 


5344 


± 


65 


0.971 


± 0.039 


162020 


4578 




68 


0, 


.800 


± 


0. 


.056 


004208 


5586 


± 


66 


0, 


.896 


± 


0. 


,044 


073526 


5615 


± 


68 


1.507 


± 0.099 


168443 


5491 


± 


69 


1, 


.610 




0. 


.042 


006434 


5741 


± 


70 


1, 


.077 


± 


0. 


,045 


074156 


5910 


± 


70 


1.667 


± 0.076 


168746 


5468 


± 


69 


1, 


.148 


± 


0. 


.049 


008574 


5942 


± 


57 


1, 


.409 


± 


0. 


.042 


075289 


6098 


± 


73 


1.245 


± 0.031 


169830 


6227 


± 


74 


1, 


.829 


± 


0. 


.042 


009826 


6184 


± 


75 


1, 


.613 


± 


0. 


.028 


075732 


5247 


± 


69 


0.934 


± 0.030 


177830 


4701 


± 


67 


3. 


.354 




0. 


.054 


010697 


5510 


± 


60 


1, 


.839 


± 


0. 


.036 


076700 


5645 


± 


69 


1.370 


± 0.047 


178911B 


5366 


± 


69 


1. 


.195 




0. 


.233 


012661 


5473 


± 


68 


1, 


.198 


± 


0. 


.041 


080606 


5461 


± 


70 


0.975 


± 0.338 


179949 


6169 


± 


73 


1, 


.194 


± 


0. 


.034 


013445 


5128 


± 


56 


0, 


.797 


± 


0. 


,025 


082943 


5952 


± 


71 


1.147 


± 0.035 


186427 


5633 


± 


71 


1, 


.190 


± 


0. 


.029 


016141 


5679 


± 


69 


1, 


.460 


d= 


0. 


,056 


083443 


5386 


± 


67 


1.055 


± 0.048 


187123 


5665 


± 


70 


1, 


.216 


± 


0, 


.043 


017051 


6269 


± 


66 


1, 


.085 


± 


0. 


,025 


089744 


6106 


± 


71 


2.182 


± 0.037 


190228 


5081 


± 


65 


2, 


.648 


± 


0, 


.057 


019994 


5999 


± 


64 


1, 


.802 


± 


0. 


.029 


092788 


5590 


± 


67 


1.075 


± 0.041 


190360 


5552 


± 


68 


1, 


.149 


± 


0. 


.028 


020367 
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Fig. 3. — (a) Teff vs (V — K) for the planet- hosting stars from the IRFM (filled circles), Ribas 
et al. (2003) (open circles) and Santos et al. (2004) (open squares). Here, K is from 2MASS 
and only stars with reliable 2MASS photometry are plotted. The solid line corresponds to 
Eq. (2). (b) Tcs vs {V — K) for the colors derived from Kurucz models for [Fe/H] = +0.5 
(dotted line), [Fe/H] = (solid line), and [Fe/H] = —0.5 (dashed line); the filled circles 
correspond to Eq. (2). 
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Fig. 4. — Comparison between the temperatures obtained in this work and those given 
by: (a) Santos et ah (2004), (b) Ribas et al. (2003), (c) Melendez & Ramirez (2003) and 
Blackwell & Lynas-Gray (1998) IRFM calibrations for the Geneva colors, (d) Sadakane et al. 
(2002) and Takeda et al. (2002). 



